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Abstract 



We have studied the color-suppressed hadronic decays of neutral B mesons 
into the final states D'^*-**'7r''. Using 9.67 million BB pairs collected with 
the CLEO detector, we observe the decays b'^^D'^tt*' and B''^D*^7r'' 
with the branching fractions ^(B^^D°7r°) = (2.74t|]J^±0.55)x 10"^ and 
B(B°-^D*07r°) = (2.20t[jJ^±0.79)xlO-^. The first error is statistical and 
the second systematic. The statistical significance of the D'^vr'^ signal is 12.1a 
{5.9a for D*'^7r'^). Utilizing the b''— >D*^*)'^7r'^ branching fractions we deter- 
mine the strong phases between isospin 1/2 and 3/2 amplitudes in the 
Dvr and D*7r final states to be cos (5/^i3=0.89ib0.08 and cos (5/^i3*=0.89ib0.08, 
respectively. 
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The decay B — >D'^*)°7r° proceeds predominantly through the internal spectator diagram 
shown in Fig. |l[ This diagram is color-suppressed, since the color of the quark-pair originat- 
ing from the W decay must match the color of the other quark pair. The rate of B —*D^*^^tt^ 
relative to 'B~^D^*^^tc~ is suppressed, crudely, by one factor of (1/3)^, and an additional 
(1/a/2)^ due to the projection of the dd state onto the pion wave function This gives a 
total suppression of 1/18 compared to the color-favored decay modes. Detailed theoretical 
calculations predict an even larger suppression of about a factor 1/50. 

So far the only established color-suppressed decays are two-body B decays into Charmo- 
nium plus neutral hadrons. A measurement of B — *-D^*^°7r° is therefore a benchmark test for 
theoretical models of hadronic B decays |jl|J^. An investigation of color-suppressed decays 
into a D meson and light neutral mesons other than a 7r° is currently underway and will be 
addressed in a future publication. 

The observation of B — >D''*'*°7r° completes the measurement of D^*^tc final states and 
allows us to extract the strong phase difference between isospin 1/2 and 3/2 amplitudes 

ii- 

In this Letter we present the observation of B —^D''*^^tc^, superseding the limits from 
our previous publication |Q. Our new analysis has significantly increased statistics and is 
based on a data sample with improved calibration and track reconstruction. The background 
shapes and the signal separation power have also been significantly improved. Charge conju- 
gates are implied throughout the paper. Our analysis uses e"'"e~ annihilation data recorded 
with the CLEO detector at the Cornell Electron Storage Ring. The integrated luminosity of 
our data sample is 9.15 fb^^ for data collected on the T(4S) (on- resonance) , corresponding to 
9.67 million BB pairs, and 4.35 fb^^ collected 60 MeV below the BB threshold (off- resonance , 
which is used for background studies. 

CLEO is a general purpose solenoidal magnet detector. Data were recorded with two 
detector configurations, CLEO II and CLEO II. V 0J§]. Cylindrical drift chambers in a 1.5T 
solenoidal magnetic field measure momentum and specific ionization (dE/dx) of charged par- 
ticles. Photons are detected using a CsI(Tl) crystal electromagnetic calorimeter, consisting 
of a barrel-shaped central part of 6144 crystals and 1656 crystals in the forward regions of 
the detector (endcaps). In the II. V configuration the innermost chamber was replaced by a 
three-layer, double-sided silicon microvertex detector, and the main drift chamber gas was 
changed from argon-ethane to a helium-propane mixture. As a result of these modifications, 
the CLEO II. V part of the data (2/3 of the total sample) has improved momentum resolution 
and particle identification. 
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FIG. 1. Diagram for the color-suppressed decay B^^D^*^^tt^. 



4 



B mesons were reconstructed by selecting high-momentum D*^*)" and 7r° mesons. Track 
quahty requirements are imposed on charged tracks and the purity of pion and kaons is im- 
proved by using dE/dx information whenever available. The 7r° candidates are reconstructed 
from isolated electromagnetic clusters of at least 30 MeV in the barrel region and 50 MeV in 
the endcaps. The mass resolution is 8 MeV in the barrel region and 10 MeV in the endcaps. 
We require that the candidate's mass is within 2.5 standard deviations (a) of the nominal 
TT^ mass. Prompt vr^'s from B decays are required to have a momentum larger than 1.8 GeV. 

D° mesons are selected in the decay modes D^^K'n'^, D°— i>K~'7r"'"7r'' and 
D^— s>K~7r"'''7r"'"7r~. The invariant mass of the D daughter particles is required to be within 2.5 
a of the known D° mass. The mass resolution depends on the decay mode and is between 
6 and 12 MeV. The momentum of the D" is required to be larger than 1.65 GeV. In the 
D'^^K~7r+7r'' mode we suppress combinatorial background by using only certain regions of 
the Dalitz plane. 

D**^ mesons are selected in the decay modes D*'^^D'^7r'' and D*''^D'^7. To reduce the 
combinatorial background in the decay mode D**^-^D°7 we require that the D" decays into 
K^7r+. We require that the mass difference m^.o-moo is within 2.5a of the known value and 
that the momentum of the D*° is larger than 1.8 GeV. The kinematic resolution of vr" and 
D^*^^ candidates is improved by a mass-constrained kinematic fit. 

B decay candidates are selected from 7r° and D*^*)° pairings that have no electromag- 
netic clusters in common. B candidates are identified using a beam-constrained mass 
Mb = yEbeam^ " Pb) whcrc Ebeam IS the beam energy and pb is the B candidate momentum, 
and an energy difference AE = Ed + Ej^o — Eboam, where Ed and E^o are the energies of the 
D(*)o and vr". The resolution in Mb depends on the D decay mode and is between 3.5 and 4 
MeV. The resolution is dominated by the beam energy spread and the tt" energy resolution. 
The resolution in AE is between 35 and 40 MeV. The energy resolution is slightly asymmetric 
due to the energy loss out of the back of the Csl crystals. The mass-constrained kinematic fit 
to the pion 4-momentum compensates for most of this effect. We accept B candidates with 
Mb above 5.24 GeV and |AE| < 300 MeV. To better suppress background from e"'"e~ — > qq 
events (continuum background), several event shape variables are combined into a Fisher 
discriminant JFd 0. For BB events (continuum background), the JFd distribution is almost 
a Gaussian and has its maximum at 0.42 (0.57). The standard deviation is 0.11 (0.12). The 
separation between the BB and continuum distributions is 1.3 a. We reject clear continuum 
events by requiring JFd< 1. For each candidate we calculate the sphericity vectors [|] of the 
B daughter particles and of the rest of the event. We require the cosine of the angle between 
these two vectors to be within -0.8 and 0.8. The distribution of this angle is strongly peaked 
at ±1 for continuum background and is nearly flat for BB events. 

The number of signal events in the sample is obtained from unbinned, extended maximum 
likelihood fits. The free parameters of the fits are the number of signal events, background 
from B decays (BB) and from continuum e"'"e~ annihilation (continuum). Four variables are 
used as input to the maximum likelihood fit: the beam-constrained mass Mb, the energy dif- 
ference AE, the Fisher Discriminant JFd, and the cosine of the decay angle of the B cos6'BHei., 
defined as the angle between the D^*-* momentum and the B flight direction calculated in the 
B rest frame. 

In each of the fits, the likelihood of the B candidate is the sum of probabilities for the 
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signal and two background hypotheses with relative weights maximizing the likelihood. The 
probability of a particular hypothesis is the product of probability density functions (PDFs) 
for each of the input variables. The PDFs for Mb are represented by a bifurcated Gaussian 



for signal, an empirical shape, Mb y (1 — x'^) exp(— — a;^)), with x = MB/Ebeam, 
for continuum and a Gaussian on top of the empirical background shapeQ for BB; the PDFs 
for AE are the sum of two Gaussians with a common mean (signal), Ist-order polynomial 
(continuum) and a sum of two Gaussians plus a Ist-order polynomial (BB); the PDFs for 
JFd are the sum of two Gaussians with a common mean; and the PDFs for cos ^BHcl. '^^^ 
2nd-order polynomials. The PDF parameters are determined from off-resonance CLEO data 
(continuum) and from high-statistics Monte-Carlo (MC) samples (signal and BB). 

Monte Carlo experiments are generated to test the fitting procedure and to obtain the 
relation between fit yield and signal branching fractions. The experiments are repeated 
several hundred times with different Monte Carlo test samples randomly selected from high- 
statistics MC samples. 

We summarize the results of the fits to CLEO data in Table |. We give results for all B 
decay modes, corresponding D decay modes and the combination of all D decay modes. We 
combine the results for different D decay modes by adding the log likelihood as a function 
of the branching fraction. Branching fractions for each mode are obtained via 

B(B«^D(*)°vr°) = ^ 

e X i3(DW) X N(BO,B^) 

The number of B° plus B , N(B°, B )=9.67 M ± 0.10 M, is derived assuming equal branching 
fractions for charged and neutral B meson decays [|10|. The uncertainty in the branching 
fractions of the T(4S) is taken into account in the systematics. The significances of the 
observed signals in the seven fits is determined from the change in -21og£ when refit with 
the signal yield constrained to zero: significance = ^2{\ogC — log £Nsig=o)- We obtain a 

total significance of 12.1cr for B^— s>D'^7r° and 5.9a for B°— s>D*°7r'^. Varying the PDF shapes 
within the systematic errors to obtain the lowest signal yield, the statistical significance 
is reduced to 9.4cr (B — >-D°7r°) and 4.2o" (B ^D*°7r°). We obtain branching fractions of 
S(B°-^DV) = (2.74;°:|±0.55)xl0-4 and B(^^D*^7r^) = (2.20t[i:| ± 0.79)xl0-^ The 
first error is statistical and the second error systematic. Our result for B^— ^D'^vr*' is higher 
than the previous CLEO upper limit We ascribe this disagreement, which is of the order 
of 3.1 cr, partly to a statistical fluctuation and partly to the description of the AE-background 
in the old CLEO publication. 

We consider sources of systematic uncertainties from the PDF shapes, D and T(4S) 
branching ratios, luminosity, possible fit bias, B candidate reconstruction and cross-feed 
between different modes. The dominant systematic uncertainty comes from the PDF shapes. 
The systematic uncertainty on the shapes is derived by varying the PDF shapes within the 
statistical errors of the fit parametrization as well as comparing the CLEO data in the AE and 
Mb sideband regions to the PDF shapes and taking differences as systematic errors. Figures 



'^The parameters of the empirical shapes for BB and Continuum background are different. 
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and show our results for B ^D^*^^7r^ with the number of signal, BB and Continuum 
background as free parameters of the fit. The fit result is projected into a signal region, 
defined in the Mb-AE plane as -0.05 < AE < 0.05 GeV, 5.275 < Mb < 5.285 GeV. The fit 
results describe the data well. The background in the sidebands is also well modeled by the 
fit. 





Mode 


Fit Yield 
(Events) 


Signifi- 
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TABLE I. Fit yields for all decay modes. Our results are based on the D branching ratios given 



in column 5 |11|. Our measurement of the B branching ratios is given in the last column. 



The observation of '^^D^*^^7t^ completes the measurement of D*^*^7r final states. This 
allows us to calculate the relative phase between the isospin 1/2 and 3/2 amplitudes in the 
D'^*)7r system. The basic relation can be expressed in an amplitude triangle: ^(D^tt"*") = 
^(D-7r+) + V2 ^(D°7r°), following the formulation in g. With the PDG values 111 
B(d\+) = (53±5)xl0-^ i3(D-7r+) = (30±4)xl0-^ S(D*V) = (46±4)xl0-^ B{D*-7i+) 
= (27.6±2.1)xl0-^ t{B+)/t{B^) = 1.073 ± 0.027, and our measurement of B°^DW07r°, 
we determine the relative phase between the isospin amplitudes to be cos^/^d =0.89±0.08 
for the Dtt final state and cos 5/_£)*=0.89±0.08 for D*7r. The ratios of isospin amplitudes 
A1/2/A3/2 are 0.70±0.11 (Dvr) and 0.74±0.08 (D*7r). A similar calculation has been per- 
formed in [0 using our preliminary results |13| and preliminary results obtained by the 



Belle collaboration [113 



Models of hadronic B decay have successfully described experimental results using two 
phenomenological parameters, ai and 02, that characterize non-factorizable contributions. 
Both are believed to be process-dependent but so far experimental data have been consis- 
tent with universal values for ai and 02. Recent work by Beneke, Buchalla, Neubert and 



Sachrajda has shown that ai is only slightly process-dependent. Based on our B^—*D^tt^ 



measurement, we derive a value a2=0.57±0.06. Comparing our result to the 02 value from 
two-body B decays to charmonium, 02=0.29 0, the process dependence of 02 is favored [|12|. 

To summarize we observed the color-suppressed decays B*^— >^D°7r° and B^— >D*°7r°. The 
number of signal events in our data sample was obtained from an unbinned extended max- 
imum likelihood fit in four variables. The measurements of the two branching fractions are 
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B(B"^DV) = (2.74l°:i±0.55)xl0-^ and -B(B"^D*V) = (2.20lg:i±0.79)xl0-^ The 
first error is statistical and the second systematic. The statistical significance of the D'^tt'^ 
signal is 12.1(7 {5.9a for D*^tt^). 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Research Corporation and the Texas Advanced 
Research Program. 
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FIG. 2. Distribution of fit input variables for B vr . The results of the unbinned, extended 
maximum likelihood fit are shown as the full line. The dotted line represents the fitted continuum 
and the dashed line is the fit result for the sum of BB and continuum background. To enhance the 
signal for display purposes, the fit results are projected into the Mb-AE signal region —0.05 < AE 
< 0.05 GeV, 5.275 < Mb < 5.285 GeV. 
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FIG. 3. Distribution of fit input variables for ^D*^'ir^ . The results of the unbinned, extended 
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